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Abstract: A general approach to accounting for retardation effects in the long-range (confining) part 
of the quark-antiquark potential is presented. The charmonium and bottomonium mass spectra are 
calculated with the systematic account of relativistic and retardation effects and the one-loop radiative 
corrections. A good fit to available experimental data on the mass spectra is obtained. 
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The relativistic properties of the quark-anti- 
quark interaction potential play an important role 
in analysing different static and dynamical char- 
acteristics of heavy mesons. The Lorentz-structu- 
re of the confining quark-antiquark interaction is 
of particular interest. In the literature there is no 
consent on this item. For a long time the scalar 
confining kernel has been considered to be the 
most appropriate one The main argument 
in favour of this choice is based on the nature of 
the heavy quark spin-orbit potential. The scalar 
potential gives a vanishing long-range magnetic 
interaction, which is in agreement with the flux 
tube picture of quark confinement of ||] , and al- 
lows to get the fine structure for heavy quarkonia 
in accord with experimental data. However, the 
calculations of electroweak decay rates of heavy 
mesons with a scalar confining potential alone 
yield results which arc in worse agreement with 
data than with a vector potential ||. The ra- 
diative Ml-transitions in quarkonia such as e. g. 
J/ip — > rjc^ are the most sensitive to the Lorentz- 
structure of the confining potential. The rela- 
tivistic corrections for these decays arising from 
vector and scalar potentials have different signs 
||. In particular, as it has been shown in ref. ||, 
agreement with experiments for these decays can 
be achieved only for a specific mixture of vec- 
tor and scalar potentials. In this context, it is 
worth noting, that the recent study of the qq 



interaction in the Wilson loop approach Q in- 
dicates that it cannot be considered as purely 
scalar. Moreover, the found structure of spin- 
independent relativistic corrections is not com- 
patible with a scalar potential. A similar conclu- 
sion has been obtained in ref. |jj on the basis of a 
Foldy-Wouthuysen reduction of the full Coulomb 
gauge Hamiltonian of QCD. There, the Lorentz- 
structure of the confinement has been found to 
be of vector nature. The scalar nature of spin 
splittings in heavy quarkonia in this approach 
is dynamically generated through the interaction 
with collective gluonic degrees of freedom. Thus 
we see that while the spin-dependent structure of 
qq interaction is well established now, the spin- 
independent part is still controversial in the liter- 
ature. The uncertainty in the Lorentz-structure 
of the confining interaction complicates the ac- 
count of retardation corrections since the rela- 
tivistic reconstruction of the static confining po- 
tential is not unique. Here we present the gen- 
eralized prescription of such reconstruction and 
discuss its implications for the heavy quarkonium 
mass spectra. 

In our preceding papers we have developed 
the relativistic quark model based on the quasipo- 
tential approach. A meson is described by the 
wave function of the bound quark-antiquark state, 
which satisfies the quasipotential equation || of 
the Schrodinger type 
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where the relativistic reduced mass is 
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E„ + 



and E a , E b are given by 
M 2 — ml + ml 



E, 
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Here M = E a + E b is the meson mass, m a , b are 
the masses of light and heavy quarks, and p is 
their relative momentum. In the centre of mass 
system the relative momentum squared on mass 
shell reads 

5 2 (M) = [M 2 - (m a + m b ) 2 ][M 2 - (m a - m b ) 2 } 



AM 2 

(4) 

The kernel V(p, q; M) in Eq. (@) is the quasi- 
potential operator of the quark-antiquark inter- 
action. It is constructed with the help of the off- 
mass-shell scattering amplitude, projected onto 
the positive energy states. Constructing the quasi- 
potential of the quark-antiquark interaction we 
have assumed that the effective interaction is the 
sum of the usual one-gluon exchange term with 
the mixture of long-range vector and scalar lin- 
ear confining potentials, where the vector confin- 
ing potential contains the Pauli interaction. The 
quasipotential is then defined by [H 



V(p,q;M) = u a (p)u b {-p)l -a.D^Q^^ 



+V v (k)T»r b;fl + V s (k)ju a (q)u b (~q), (5) 

where as is the QCD coupling constant, is 
the gluon propagator in the Coulomb gauge and 
k = p — q; 7^ and u{p) are the Dirac matrices 
and spinors with e(p) = \/p 2 + m 2 ■ The effective 
long-range vector vertex is given by 



r M (k) = 7^ + l^ a ^k v 



(6) 



where k is the Pauli interaction constant charac- 
terizing the anomalous chromomagnetic moment 



of quarks. Vector and scalar confining potentials 
in the nonrelativistic limit reduce to 

V v (r) = (1 - e)Ar + B, V s (r) = sAr, (7) 

reproducing 

Konf(r) = V s (r) + V v (r) = Ar + B, (8) 

where e is the mixing coefficient. 

The retardation contribution to the one-gluon 
exchange part of the qq potential is well known. 
For the confining part of the qq potential the re- 
tardation contribution is much more indefinite. 
It is a consequence of our poor knowledge of the 
confining potential especially in what concerns 
its relativistic properties: the Lorentz structure 
(scalar, vector, etc.) and the dependence on the 
covariant variables such as k 2 = k 2 , — k 2 . Never- 
theless we can perform some general considera- 
tions and then apply them to a particular case of 
the linearly rising potential. To this end we note 
that for any nonrelativistic potential V(— k 2 ) the 
simplest relativistic generalization is to replace it 
by F(fc 2 -k 2 ). 

In the case of the Lorentz-vector confining 
potential we can use the same approach as for the 
one-gluon exchange even with more general ver- 
tices containing the Pauli terms, since the mass- 
shell vector currents are conserved here as well. 
It is possible to introduce alongside with the "di- 
agonal gauge" the so-called "instantaneous gauge" 
|J which is the generalization of the Coulomb 
gauge: 

V v {kl - k 2 )u a (p)u b (-p)r^r bM M a (q)M b (-q) 

= z2 a (p)w b (- P )|^(-k 2 )r^ 
-[v v (~k 2 )r a ■ r b + ^(-k 2 )(r a • k) 

x(r b -k)] J Ua (q) Wb (-q), (9) 
where 

V v (k 2 - k 2 ) = V v (-k 2 ) + k 2 V{s(-k 2 ) 

and 

k 2 = ( eo (p)-e (q))( eb (q)- eb (p)) - 

°(10) 

with the correct Dirac limit in which the retarda- 
tion contribution vanishes when one of the par- 
ticles becomes infinitely heavy |fOf. 
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For the case of the Lorentz-scalar potential 
we can make the same expansion in k 2 , which 
yields 



Vs{kl 



V s (-k 2 ) + k 2 V!,(-k 2 



(11) 



But in this case we have no reasons to fix k 2 in 



the only way (10). The other possibility is to 
take a half sum instead of a symmetrized prod- 
uct, namely to set (see e. g. fll], 



kl = -[Mp)-£ a (q)) 2 + Mq)-e h (p)) 2 



= > 2 -q) 5 



(12) 



The Dirac limit is not fulfilled by this choice, but 
this cannot serve as a decisive argument. Thus 
the most general expression for the energy trans- 
fer squared, which incorporates both possibilities 
@ and @ has the form 

kl = A(e„(p) - e a (q))(€ fe (q) - e 6 (p)) + (1 - A) 
x^[(ea(p)- e a(q)) 2 + ( eb (q)- £6 (p)) 2 ] 

a A (p 2 -q 2 ) 2 

Am a m b 

+(l-^(P 2 -<0 2 (4 + ^), (13) 

where A is the mixing parameter. 

Thus the spin-independent part of qq poten- 
tial with the account of retardation corrections 
takes the form: 

Vsi{r) = V c (r) + V conf (r) + V VD {r) 

' 1 V + A[Vfc(r) + (l + 2 K )W(r)],(14) 



where the velocity-dependent part 

Vy D (r) = V&(r) + V&(r) + V&(r), 



(15) 
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and {. . .}w denotes the Weyl ordering of opera- 
tors. Making the natural decomposition 



VVo(r) 



1 



m a m b 
mi mi; 



p 2 V bc (r) + ^p-V c (r) 



(P • if 



w 
w 



we obtain [Q for the corresponding structures 
with Ay = 1 and including one-loop radiative 
corrections in MS renormalization scheme: 



Vc(r) = -- 


V bc (r 










+ 




Vc(r) 
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ln(/xr) 



2vr 

„2( ,,T\ 



r 



3 2tt 
Ar + B, 

2 /3 a 2 (M 2 ) 

3 2tt 
1 — e eXs 

2 + ~Y 



Ar, 



V de (r) = --(l + \ s )Ar + B, 
V e (r) = -~(l-\ s )Ar, 



(17) 



where 

a v (n 2 ) = ctsip 2 ) 
31 10 



1 



1ePo\ ai s {n 2 ) 



ai = rif, 

3 9 ; ' 



Oi 

4 



A) = 11 - g"/- 



Here n./ is a number of flavours and is a renor- 
malization scale. It is easy to check that the ex- 
act Barchielli, Brambilla, Prosperi relations p3| 
following from the Lorentz invariance of the Wil- 
son loop 

V de -^V bc + ^(V c + V o ) = 0, 

K+ ^ + r_m^^ (18) 

are exactly satisfied. 

The expression for spin-dependent part of 
the quark-antiquark potential with the inclusion 
of radiative corrections can be found in rcf. jl2| . 
Now we can calculate the mass spectra of heavy 
quarkonia with the account of all relativistic cor- 
rections (including retardation effects) of order 
v 2 jc 2 and one-loop radiative corrections. For this 
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purpose we substitute the quasipotential which is 
a sum of the spin-independent and spin-dependent 
parts into the quasipotential equation. Then we 
multiply the resulting expression from the left 
by the quasipotential wave function of a bound 
state and integrate with respect to the relative 
momentum. Taking into account the accuracy of 
the calculations, we can use for the resulting ma- 
trix elements the wave functions of Eq. (|l|) with 
the static potential 



V NR (r) 



4 ay(/i 2 ) 
3 r 



Ar + B. 



(19) 



As a result we obtain the mass formula [m a 
rrib = m) 



b 2 {M) 



+ A(S -r)(S b T 



W +(a)(L-S) + (b)( 



-(S„-S 6 ) 



} + (c)(S Q -S b ), (20) 



where the first term on the right-hand side of the 
mass formula contains all spin-independent con- 
tributions, the second term describes the spin- 
orbit interaction, the third term is responsible for 
the tensor interaction, while the last term gives 
the spin-spin interaction. 

To proceed further we need to discuss the 
parameters of our model. There is the follow- 
ing set of parameters: the quark masses (mj and 
m c ), the QCD constant A and renormalization 
point /i in the short-range part of the QQ po- 
tential, the slope A and intercept B of the linear 
confining potential (||), the mixing coefficient e 
(0), the Ion g-range anomalous chromomagnetic 
moment k of the quark (||), and the mixing pa- 
rameter As in the retardation correction for the 
scalar confining potential. We can fix the values 
of the parameters e = — 1 and k = — 1 from the 
consideration of radiative decays || and compar- 
ison of the heavy quark expansion in our model 
|l4| , 15 with the predictions of the heavy quark 
effective theory. We fix the slope of the linear 
confining potential A = 0.18 GeV 2 which is a 
rather adopted value. In order to reduce the 
number of independent parameters we assume 
that the renormalization scale /i in the strong 
coupling constant a s (/i 2 ) is equal to the quark 
mass. We also varied the quark masses in a rea- 
sonable range for the constituent quark masses. 
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Table 1: Charmonium mass spectrum. 

The numerical analysis and comparison with ex- 
perimental data lead to the following values of 
our model parameters: m c = 1.55 GeV, nib = 
4.88 GeV, A = 0.18 GeV 2 , B = -0.16 GeV, fi = 
m Q , A = 0.178 GeV, e = -1, k = -1, X s = 0. 
The quark masses m c ^ have usual values for con- 
stituent quark models and coincide with those 
chosen in our previous analysis [||. The above 
value of the retardation parameter As for the 
scalar confining potential coincides with the min- 
imal area low and flux tube models jl6| , with lat- 
tice results jffl and Gromes suggestion jll| . The 
found value for the QCD parameter A gives the 
following values for the strong coupling constants 
a s (m 2 c ) w 0.32 and a s (m 2 b ) w 0.22. 

The results of our numerical calculations of 
the mass spectra of charmonium and bottomo- 
nium (in GeV) are presented in Tables [j] and ||. 
We see that the calculated masses agree with ex- 
perimental values within few MeV and this dif- 
ference is compatible with the estimates of the 
higher order corrections in v 2 /c 2 and a s . The 
model reproduces correctly both the positions 
of the centres of gravity of the levels and their 
fine and hyperfine splitting. Note that the good 
agreement of the calculated mass spectra with 
experimental data is achieved by systematic ac- 
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State Particle Theory Experiment [lq 



1% 

l 3 ^ 
1 3 P 

l 3 Fi 
1 3 P 2 
2% 
2 3 S 1 

1 3 L>! 
l 3 D 2 

1 3 D, 



'lb 
T 
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9.400 
9.460 
9.864 
9.892 
9.912 
9.990 
10.020 
10.151 
10.157 
10.160 



9.46037 
9.8598 
9.8919 
9.9132 

10.023 



influence the wave functions. As a result the de- 
cay matrix elements involving heavy quarkonium 
states remain mostly unchanged. 

We are grateful to the organizers for the nice 
meeting and stimulating discussions. Two of us 
(R.N.F and V.O.G.) were supported in part by 
Russian Foundation for Fundamental Research 
under Grant No. 00-02-17768. 
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